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TECHNICAL PAPER 


DAMPING SEALS FOR TURBOMACHINERY 
INTRODUCTION 


The stable operation of turbomachinery depends mainly on the damping of the rotor motion [ 1 ] . 
Squeeze film dampers are designed for this very purpose [ 2 ], and also shaft seals are employed [3,4] . Seals 
are most convenient, since they are already included in the turbomachine design. However, the damping of 
a seal also generates a whirl driver force as shown in the pioneering work of Henry Black [5] who died 
January 17, 1980. He and his service to all are well-remembered [6] . 

The seal damping, normally effective only to twice the first critical speed, becomes an unstable 
whirl driver at higher speeds. Thus, high performance turbomachinery which has to run above the first 
critical speed is limited by an upper boiind twice the first critical speed unless additional damping is applied. 
The fluid shear forces between the rotating journal and the stationary bushing produce a circumferential 
bulk flow, known as Couette flow [7], at half the shaft speed. The fundamental relation between Couette 
flow and whirl forces was recognized by Black [5]. Black, Allaire, and Barrett [8] found for the circum- 
ferential flow, as an approximation, an exponential function with the axial position as the independent 
variable, the inlet swirl as the boundary condition, and the Couette flow at half speed as the asymptote. 
It is shown that the whirl force can be reduced with a low or reversed inlet swirl under the assumption that 
the swirl can be transferred from the inlet cavity to the narrow seal gap. 

This paper proposes the application of the seal surface roughness for stabilization and sealing. The 
distribution of the circumferential flow is approximated by the asymptotic state of the Couette flow which 
is assumed to prevail over the total seal length. The approximation is justified by the results of Black et al. 
[8], who show that seal friction advances the asymptotic approach upstream toward the inlet. The formu- 
lation treats the Couette flow generation by surfaces of different roughness without invoking a dominant 
axial flow as in References 5 and 8. This generalization is more representative of long seals. The turbulent 
seal friction is given by Moody’s pipe friction formula [9] for Reynolds numbers up to 10^ and above. 
Linearization is applied at the steady flow point to aid the integration of the flow effects over the cylin- 
drical surfaces. The steady flow is determined with the Newton-Raphson’s iteration method. The dynamic 
parameters are derived for small deflections about the seal’s center position. Numerical results are given 
for the oxygen and the hydrogen turbopumps of the Space Shuttle main engine. 


ROTOR STABILITY 


This section briefly reviews the purpose of the rotor damping with the help of a simple model. The 
model consists of a rotor (Fig.* 1) with a mass m^, eccentricity eh^, combined bearing and seal stiffnesses Ky 

and K^., seal damping C, speed ft, and Couette flow velocity v. The whirl cross coupling is given by cftC = 

(v/r)C, with the Couette factor c and the seal radius r. The centrifugal force caused by the unbalance is 
given in equation (1). The components of the force are obtained with matrix of the equation (2) that 

rotates the unit vector i of equation (3) by the angle a. The force is balanced by the dynamic stiffness 
effect given in equation (4). The inverse is the response of equation (5). 



ROTOR WITH SPEED n, mass, 
m^, AND ECCENTRICITY, eh^ 

STATOR 



Figure 1. Simple rotor model. 
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The common denominator of equation (5), given in equation (6), results from the determinant of 
equation (4). The characteristic equation D = 0 determines the stability. The mapping of the contour of the 
right half of the complex s-plane (see Nyquist criterion [10]) yields the stability criterion. The limit is 
found by rearranging equation (6) for s = jco and setting D = 0. The imaginary part of D yields the critical 
speed and rotor resonance [Eq. (7)] ; and the real part of D, the stability limit (inequalities 8 to 10). 
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cOq = V (Ky + K^)/2m^ 

(7) 

Ky+K^ 


K K^+ (cfiC)2 + (Ky+ 14 < [Cy + m^(K^+ K^)] 

zm^ 

(8) 

Ky +K^ {Ky-K^)2 

(cnC)2 < C" ^ z ^ y z 

(9) 

2m^ 4 


SI < o^q/c 
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The stability condition of inequality (9) is obtained by rearranging the inequality (8). The special 
case of a uniform suspension (Ky = K^) yields the simplest stability criterion : inequality (10). The latter 

satisfies also the inequality (9) on the conservative side. Inequality (10) limits the speed of turbomachinery 
that operates above the critical speed. For example, a half-speed Couette flow has a c = 0.5 and limits the 
speed at twice the critical speed (J2 < 2co^) as, e.g., early in the development of the high pressure hydrogen 

turbopump of the Space Shuttle main engine [4]. If the Couette flow has a c = 0.25, then the Umit is 
extended to ^2 < 4 cOq. The objective of this study is to define a seal configuration with a Couette factor 
c<0.5. 


SEAL GEOMETRY 


Figure 2 shows the cross section of the seal for an exaggerated gap situation. The eccentricity eh^ 

(0 < e < 1) points with the time dependent angle a:(t), and the gap is measured at an angle j3. Both angles 
are inertially referenced. The cylindrical coordinates are x (axial position) and q = rjS (r = R circumferential 
position). Equations (11) through (15) give, in respective order, the rotational transformation matrix Q for 
the angle ]3-a(t), the right angle rotation matrix II, the unit vector i, the partial derivative after j3, and the 
partial derivative after time. Equations (16) through (21) are, in respective order, the compact expressions 
of cosine and sine, seal gap h, the differential relative to the gap, the circumferential partial derivative rela- 
tive to the gap, and the time partial derivative relative to the gap. Equation (18) follows from Figure 2, 
R - r = h^, R = r, and h = R - [r + ehQCos(|3 - a)] = h^[ 1 - e cos(j3 - a)] . 



AXIAL BULK FLOW VELOCITY 
RADIAL GAP h 


ROTOR WITH RADIUS r, SHEAR FORCE r , 

AND SURFACE VELOCITY w 

STATOR WITH RADIUS R AND SHEAR FORCE r. 


COUETTE FLOW VELOCITY 


Figure 2. Seal cross section. 


3 



( 11 ) 


Q = 


n = 


cos(j3-a) -sin(/3-o:) 
sin(/3-a) cos(|3-a) 

b -1 


1 0 


I = 


1 0 
0 1 


= -n^ 


1 = 


iT = [1 0] 


9Q 

— = ng 
dp 


9Q 

^ = -cd3Q 

at 


cos(j3-a) = i^Qi 


sin(/3-a) = -i"^IIQi 


h = 11 q( 1 - eilQi) 


dh 


P = -ei^Qi = -E 

1^0 


1 9h -e .T_^. 
— — =-i^UQi 
ho aq r 


1 ah T 

_ = . aelDQi 

ho at 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 


FLOW FRICTION 


The flow frictinn is given by the pressure drop in pipes Ap = p(u^ /2)(L/D)f with the friction factor 
f [9] , length L, hydraulic diameter D (D = 2h for seals with radial gap h), roughness k (k/h roughness relative 
to gap h), axial bulk flow velocity u, and fluid density p [11]. The seal has the axial and circumferential 
bulk flow velocities u and v and the rotor surface velocity w (Fig. 2). The fluid shear stresses on the surfaces 
of the stator (bushing) and the rotor (journal) are Tg in equation (22) and Tj. in equation (23). The pressure 

gradient is given in equations (24) and (25) with the axial and circumferential components. The Reynolds 
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numbers are given in equations (26) and (27), the frictions in equations (28) and (29), the friction s Taylor 
expansions in equations (30) and (31), the friction velocity factors in equations (32) and (33), and the 
friction gap factors in equations (34) and (35); all equation pairs are for the stator and rotor, respectively. 
The fluid viscosity is ju in equations (26) and (27) and the surface roughnesses kg and kj. in equations (28), 
(29), and (32) through (35) for the stator and rotor, respectively. 
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( 34 ) 


hg = ^I^^dO-^kg/h+lOVRs)^^^ 

hj. = 2:2^ (lO'^yh + (35) 


COUETTE FLOW 


The Couette flow is generated by the rotor shear force and resisted by the stator shear force. An 
eccentricity causes an additional circumferential flow that produces the pressure gradient component of 
equation (25). The latter is assumed to be negligible, as expressed by the equation (36), which yields the 
Couette factor c = v/w. The limits given in equations (37) and (38) follow from the circumferential and the 
axial flow dominances. The Couette factor varies little with u, as Figure 3 shows. 


■gN/u^+v^ + (v-w) fj.N/u^+(V“W)^ =0 

(36) 

= 0; c = — = 1/(1 +s/W 
w ^ 

(37) 

= OO ; c = — = 1/(1 + fg/fj.) 

(38) 


w 


w/w 0.5 



Figure 3. Couette factor c = v/w versus the axial flow velocity 
ratio u/w for friction ratios fg/fj.- 
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FLOW CONTINUITY 


The motion of the gap pumps the assumed incompressible fluid, according to equation (39), to pre- 
serve the volume. The equation is further simplified to equation (40) by canceling terms, assuming no gap 
taper (9h/3x = 0), and by assuming a constant circumferential flow (0v/8q = 0). Equation (41) is obtained 
by substituting equations (20) and (21) in equation (40), and by approximating 1/h with l/hg. The intro- 
duction of the definitions of the equations (42) and (43) and the integration of equation (41) yield the axial 
velocity of equation (44) with uj as the inlet velocity. 


3u 9v 9h 9h 

uhdq + vhdx = uhdq + h — dxdq + vhdx + h — dqdx + v — dqdx + — dxdq 

9x 9q 9q 9t 


(39) 


9u ^ v9h Q 

9x h9q h9t 

9u T . . V X 

— = i^(el - aeEDQi + - ei^OQi = A 
9x r 

y = q ;._ = a-Cii 

r 

A = i^(el - 7eEl)Qi 
u = Uj + xA 


(40) 

(41) 

(42) 

(43) 

(44) 


PRESSURE DISTRIBUTION 


The dynamic pressure gradient of equation (45) follows from equation (24) by adding the mass 
acceleration of equation (46). The equation (47) is obtained by expanding all factors of equation (45) with 
first order terms and by utilizing the expanded friction factors of equations (30) and (31). The multiplica- 
tion of all factors and the neglecting of second order terms yield equation (48), which is the Taylor expan- 
sion of equation (45). The definitions of equations (49) through (51) lead to the compact form of the 
Taylor expansion, equation (52). The time derivative of equation (46) applied to equation (44) yields equa- 
tion (53). The definitions of equations (42) and (56) result in a further simplification, equation (54). The 
pressure distribution of equation (58) is obtained by integrating equation (52); by substituting the defini- 
tions of equations (19), (43), (55), and (56); and by employing the approximation of equation (57). 


17 " ^ [fsx/u^+v^ +fj.yu2+(v-w)2] + p^ (45) 
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( 19 ) 



Auj = u - uj = xA 
B = i’^[(e-72e)I-(ae+27e)njQi 

-ui^ 

— = ujCu-uj) = UJXA 

x2 du j 

pj - p = xujF + — A(F+G) + xujE(F+H) + xp +— pB + xpujA 


(55) 

(56) 

(57) 

(58) 


BOUNDARY CONDITIONS 

The pressure drop over the total seal gap length L of equation (59) follows from x = Lin equation 
(58). The inlet loss fj reduces the inlet cavity pressure p^ to the inlet pressure pj of equation (60), assuming 

Uq = 0. The outlet pressure P 2 equals the outlet cavity pressure P 3 of equation (61), assuming a sudden 
discharge into a large outlet cavity (f 2 = 1) and U 3 = 0. Equation (62) gives the pressure drop from one to 
the other cavity. Equation (63) describes the inlet flow velocity u^ as the sum of the steady (u^) and the 
oscillatory (uj^) velocities. Equation (64) results from equations (62) and (63), u^ » U|j, ujF = u^F + 
Ujj(F+G), and from neglecting second order terms. The steady-state part of equation (64) is given in equation 

(65) as a Taylor expansion, and in equation (66) in the iterative form after Newton-Raphson. Equation (67) 
is the oscillatory component of equation (64) with s replacing d(-)/dt and tacitly assuming the Laplace 
transform for the variables. Equation (68) defines the time constant r, and equation (69) the limit Tq that is 

found by conservatively neglecting pu 2 (l+fj). Equation (70) follows from equations (58) and (60). Equa- 
tion (71) follows from equations (63) and (70), u^ » u^, U 2^/2 - u^^/2 u^u^^, ujF = u^F -t- u^,(F+G), 

and by neglecting second order terms. The pressure distribution of equation (71), together with equations 

( 66 ) and (67), expresses the effects of the boundary conditions and the seal gap. 

L 2 du 2 j ^2 

Pi - p 2 = LujF -I- LujE(F+H) + — A(F+G) + Lp + -^ pB + LpujA (59) 
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Ui 


Pq-P = Pi + 


= p E(F+H) + pA + sp] + [A(F+G) + pB] 


P = Pq-P— (1+fi) - xUgF - xUg [E(F+H) + pA] - y [A(F+G) + pB] 


- uj, [pUgC 1 +f 2 ) + x(F+G+sp) ] (7 1 ) 


DYNAMIC SEAL PARAMETERS 


The fluid forces are obtained by integrating the pressure distribution of equation (71). First, the 
axial integration gives the radial force in equations (72), (75), (81), and (87); and second, the circumferential 
integration gives the total fluid force in equations (88) and (107). Equation (73) defines the average pressure 
p^ which does not contribute to the side force. Equation (74) defines the factor N of the variable velocity 

Ujj as derived from equations (67) and (72). Equation (74) is further simphfied by substituting with equa- 
tions (68) and (69). The substitution of equation (67) in equation (72) yields equation (75) and the defini- 
tions of equations (76) through (84). The time constants of equations (68), (69), (80), (84), and (85) are 
ordered like the inequalities (86). The upper limit of equation (85) is derived by selecting a v/u ratio that 
minimizes the stator and the rotor frictions in F+G [equations (49), (50), and (69)]. The circumferential 
integration of equation (88) is vectorially accomplished to advantageously use the matrices of equations 
(89) and (90). The integration of the functions E, A, and B [equations (19), (43), and (56)] involves the 
function type of equation (9 1 ). The latter employs equation (89) to separate the angles a and (3 via the 
matrix product. The integrations of equations (92) lead to a diadic that, when integrated, yields a simple 
result. Equations (93) through (95) relate the polar and rectangular coordinates and their derivatives. Equa- 
tions (96) through (98), the integrals of the functions E, A, and B, are transformed with equations (42) and 
(93) through (95); and s = d(*)/dt to rectangular coordinates, tacitly assuming Laplace transformed variables. 
Equations (87) and (96) through (98) yield the definitions of equations (99) through (107). Equations 
(102) through (106) give the frequency response of the dynamic parameters. The frequency response is due 
to the time constants t, t^, and tj,, which are small as indicated by the limit of equation (85) and inequal- 
ities (86). The frequency response is obtained with s = jco and by replacing l/(l+sr) = (l-sr)/(l+cj^r^) and 
s^ = -sco^ . The (cot)^ terms are usually negligible; however, (cS2)^ m^. and tKq are not [equations (102) and 

(107)]. The vector equation (107) is the dynamic model of the seal. The whirl driver is the off-diagonal 
elements that are proportional to the Couette factor c and the rotor speed S2. 
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2(1 +st) 


(76) 


T /I 1 \ T 


(77) 


n-1/3 = - 


b 


2 1+sr 


(78) 


b = 


’’o 3 \r Tq/ 3rt 


(79) 


’’b 


pL 


4pua(l+fi) + L(F+G) 


(80) 


r . ^ F+H A 

/ pdx = Lpa + E — u„d ■ + - 

J ^ 2 ^ 1+sr 1+sr 2 

o 


pu^d + 1 (F+G) b( 1 +STb)| + B ^ pb (8 1 ) 


•j 


L3 1+STb 

— pb 

4 1+sr 


PUad + -(F+G)b(l+STb) = r 


= r f ^a' ( 1+^1 ) + (F+G) + s j (F+G) 


= a(l+STf,) 


(82) 


r F+G 
a = -.L-— 
’’a 6 


( 83 ) 
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pi? (F+G) 

’’a " pu^d+fi) [6pua + 4L(F+G) + V (F+G)^ 


( 84 ) 


4h. 


r„ = 




(85) 


’’u > '^o > “T > Tb > Ta 


( 86 ) 


V F+H \? 


/ pdx = Lp„ + E — u„d • — — + A — a + B — pb 

J ^ 2 ^ 1+sr 2 1+ST 4 1+ST 


(87) 


27t / L \ 

r f I f pdx j id|3 

o \ o / 


(88) 


Q = Q« Q|3 


(89) 


Q|3 


cos/3 -sinj3 

sinjS cos|3 


(90) 


i"*"(aQl“bQlI)Qi = cos(]8-<i(:) + sin(j3-a) = 


L^oJ 


i^QjOa 


L^o. 


(91) 


2it 2it 

J iT(aoI-boE)Qi • idjS = / ii^Q^ d/3 • Q„ 
o o 


27T 

/ 


cos^ P cos/3sin/3 
sin/3cos/3 sin^ /3 


d^Qa 


= TTQa 


= 7r(aQl + boII)Qj3^i 


(92) 


[:] 




( 93 ) 
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= hQ(eI+d:en)Q£^ 


(94) 


= 11 q [ (e-a^ e)I + (ae+2ae)ll ] 


2 tt P -I 

ho / = ’ThoCQ^i = Trj^^J 


(95) 


(96) 


2 tt r 

11q / AQ^idjS = 7rhQ(eI+7eII)QQ4 ^ ^ 

o L « 


- 7Tcr2n 


= 7r(sI-cJ2II) 


(97) 


2ir 


ho f BQp id|3 = 7rho[(e-7^ e)I + (ae+2-ye)II] Q^,! 


7T 


“ — 


” “ 


“ “ 


- - 

y 

z 

- 7r(ci2)^ 

y 

z 

- 25tcJ2I1 

y 

z 

= 7t[(s 2 -(cJ2)2)I-s2cS2H] 

y 

z 

- « 


(98) 


7rr L 


Co = ^^^(F^G) 

hoTal2 


m 


o 


7TT T 



K = 


K, 


1 H-CO^T^ 


(99) 


( 100 ) 


( 101 ) 


( 102 ) 
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(103) 


C = 


Cq(1+co^ TTg) + m^co^ - tKq 


1+OJ^T^ 


111^(1 +co^TT|3) - C^Cr-Ta) 


m = 


IWt^ 


Cc = 


CoCl+co^TTa) + (t-t^) 


l+wV^ 


(104) 


(105) 


m^Cl+co^TTb) 
" l+coV^ 


(106) 



K-(cJ2)^ + sC + m 

c^2[C(, + s(m+nij,)] 


y 

-cS2[C(, + s(m+m^)] 

K-(cJ7)^ + sC + m 


z 


(107) 


STABILITY 


The dynamic models of turbomachines are far more complex than the simple model of Figure 1. 
The dynamic models consist of a sum of normal modes with gyroscopic moments and fluid forces added. 
The normal modes are the sum of resonators. The resonators are tuned to individual resonances, such as 
the critical speeds of the rotor. At one critical speed, one resonator dominates; and thus, the simple model 
of Figure 1 appears to be representative of the rotor behavior when the rotor mass is replaced by the gen- 
eralized mass mg and the suspension stiffness by the generalized stiffness Kg at the seal location. The 

stability is assessed by adding Kg + s^mg to the diagonal of the matrix in equation (107). The combined 

parameters are defined in equations (108) and (109). The determinant of the modified matrix of equation 
(107) gives the characteristic equation (1 10). (See Eq. 6.) The imaginary part of D = 0 yields the resonance 
of equation (111); and the real part, the stability limit of the inequalities (112) through (114). The ratio 
^l/cOo of equation (115) is approximately 1/c, because C = and (m+m^^) is small. In other words, 

the Couette factor controls the speed limit. 


K^ = K + Kg (108) 

m^. = m + mg (109) 
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D = [Kj - (cS2)2 m^.] 2 + (cny + s2 [C(K^ - (c^)^ m^,) + (cn)^ C^.(ni+mc)] 


+ [C^ +2(Kj. - (cJ2)^ nig)nij. + (cJ2)* (m+rrip)* ] + 2mj.C + s'* = 0 


K m„ 

cOq2 = (cS2)2 

^ rrij. nij. 


2[Kj.-(cn)2mc]^ +(cn)2Cc^ < [Kr-(cn)2mc]C2/mr+2[Kr-(cn)2mc]2 

+ [Kj. - (cJ2)^ ni(,] (cfl)^ (m+m^)^ /rrij. 

(c£2)^C^^ < cOq^ [C^ + (cJ2)^ (m+m^,)^ ] 

C 

n < ^ 

^ \Z^c 

^ C/c 

*^o y/C^ - cOq^ (m+nij,)^ 


( 110 ) 

( 111 ) 

( 112 ) 

(113) 

(114) 

(115) 


NUMERICAL RESULTS 


Table 1 compares a Rocketdyne analysis after Black for a uniform Couette flow factor c = 0.5 with 
the new analysis for equally polished seal surfaces of relative roughnesses kg/h^ = kj./hp = 0.002. The results 

are similar and differ mostly at the fluid film stiffness. The difference is most likely due to Black’s assump- 
tion of a dominant axial flow, while the new formulation is not thus restricted. Table 2 lists the recent 
design parameters for seals of the high pressure fuel turbopump’s (HPFTP) first and second pump inter- 
stages, and the high pressure oxidizer turbopump’s (HPOTP) inducer shroud. Table 3 gives the results for 
damping seals with rough strators of kg/h^ = 0.2 and polished rotors of kj./hQ = 0.002. A high roughness 

ratio is required because of the cubic root in Moody’s friction equations (28) and (29). Table 3 demonstrates 
that speed limits of 2.8 and 3.9 times the first critical speed are feasible. The ratios are calculated 

with equation (115) for first critical speeds of 18,000 RPM and 15,000 RPM of the HPFTP and the HPOTP. 
respectively. An inlet loss factor of fj = 0.5 is applied in all cases. A comparison of Tables 1 and 3 shows 

that damping seals also have a lower leakage velocity u (70 percent). 
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TABLE 1 . COMPARISON BETWEEN THE BLACK AND THE NEW MODEL 


Dimensions 

HPFTP 

Interstage 

Seals 

Black Model 
Rocketdyne 
Data 7/13/76 

New Model 


103 Ib/in. 

Ib-s/in. 

lb 

ft/s 

1 







c 

K 

c£2Cg 

C 

cf2mj. 

m 

u 


285 

209 

107 

10.7 

2.11 

1004 

2 

404 

205 

93 

i 

6.5 

0.56 

1022 

2 


The comparison is for the straight smooth seal with a 3.14-in. diameter, 1.5 in. length, 37,360 
RPM, and a relative roughness of kg/h^ = k^/hg = 0.002 for the stator and the rotor. 


TABLE 2. DESIGN PARAMETERS OF SEALS FOR THE HIGH PRESSURE FUEL 
TURBOPUMP (HPFTP) AND THE HIGH PRESSURE OXIDIZER TURBOPUMP 
(HPOTP) OF THE SPACE SHUTTLE MAIN ENGINE (SSME) 


1^ Dimensions 


inch 


psia 

RPM 

Ib/cft 

Ib/ft-s 

Seals 

Gap 

Dia. 

Length 

In 

Out 

Speed 

Density 

Viscosity 

HPFTP First 
Interstage 

0.0055 

3.56 

1.8 

2432 

288 

36860 

4.64 

1.06 • 10-5 

HPFTP Second 
Interstage 

0.0055 

3.56 

1.8 

4810 

2557 

36860 

4.33 

0.86 • 10'5 

HPOTP Inducer 
Shroud 

0.01 

4.70 

1.5 

1824 

433 

30375 

64.5 

8.00 • 10-5 

HPOTP Inducer 
Shroud 

0.02 

4.70 

1.5 

1824 

433 

30375 

64.5 

8.00 • 10-5 
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TABLES. DAMPING SEAL PARAMETERS 


Dimensions 

lO^lb/in. 

Ib-s/in. 

lb 


ft/s 


10^ 


10'® 

s 

10'3 

1 


Seals 

K 

cnC^ 

C 

cnm^ 

m 

u 

V 

w 

Rs 

Rr 

T 

’"b 

’’a 


c 

^0 

HPFTP First 
Interstage 

340 

182 

182 

4.9 

1.3 

660 

142 

573 

271 

316 

24 

16 

15 

2 

4.0 

3.9 

HPFTP Second 
Interstage 

358 

179 

181 

4.5 

1.2 

703 

140 

573 

333 

383 

23 

15 

14 

2 

4.1 

3.9 

HPOTP Inducer 
Shroud, 0.01 in. 
Gap 

168 

240 

202 

37 

7.3 

154 

204 

623 

343 

600 

137 

78 

65 

46 

3.1 

2.9 

HPOTP Inducer 
Shroud, 0.02 in. 
Gap 

116 

123 

101 

23 

3.2 

217 

195 

623 

783 

1290 

166 

72 

46 

68 

3.2 

2.8 


The stator and rotor roughnesses are kJhQ = 0.2 and kj./hQ = 0.002, respectively, cj = 0; and the critical speeds for the 
HPFTP and the HPOTP are 18,000 RPM and 15,000 RPM, respectively. 


CONCLUSION 


The dynamic and static seal parameters of seals with rough surfaces are derived without invoking a 
dominant axial flow. Whirl stability and leakage are shown to be controllable by employing a stator with a 
high surface roughness. The feasibility of damping seals is analytically demonstrated, but experimental 
verification is needed. The proposed seals are simple and should be readily applicable to turbomachinery. 
Speed limits can thus be raised, rotor vibrations reduced, bearing life improved, and costly shutdowns 
avoided. 
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